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I 

19.  Abstract  (Continued)  j 


The  maximum  low  frequency  decoupling  factor  in  salt  is  approximately  100.  In  a  17  meter 
air-filled  salt  cavity,  the  explosion  is  fully  decoupled  at  yields  less  than  0.2  kilotons,  slightly  less  ' 
than  the  yield  predicted  using  the  Latter  criterion.  A  low  frequency  decoupling  factor  of  j 
approximately  50  is  predicted  for  a  1.5  kiloton  explosion  with  the  same  material  and  geometry. 


TAn  important  result  of  this  research  is  that  the  low  frequency  source  amplitude  increases  ; 
very  slowly  as  the  yield  is  increased  above  the  yield  required  for  full  decoupling.  In  the  1 
materials  studied,  the  decoupling  factor  decreases  by  only  about  a  factor  of  two  at  a  yield  j 
eight  times  the  Latter  decoupling  criterion.  These  results  suggest,  for  example,  that  it  would  j 
be  possible  to  detonate  a  24  KT  explosion  in  a  cavity  sufficient  to  fully  decouple  only  a  3  KT  i 
explosion,  and  still  achieve  a  decoupling  ratio  half  as  great  as  would  be  obtained  from  the  I 
much  larger  cavity  that  would  be  required  to  fully  decouple  the  larger  explosion.  j 


At  high  frequencies,  there  is  little  difference  between  simulations  with  linear  and  nonlinear 
material  models  for  decoupled  or  partially  coupled  explosions,  except  for  attenuation  due  to 
pore  crushing  All  simulations  show  a  strong  high  frequency  peak  due  to  cavity 
reverberations^  The  dominant  frequency  of  this  peak  is  a  function  of  yield,  increasing  from 
about  100  Hz  for  the  overdecoupled  simulations  to  about  500  Hz  for  partially  coupled 
simulations  for  the  yields  and  cavity  sizes  studied  here.  Decoupling  factors  decrease  as  a 
function  of  frequency,  however  this  decrease  starts  at  higher  frequencies  for  smaller, 
overdecoupled  explosions  than  for  larger,  fully  decoupled  and  partially  coupled  explosions. 
As  a  result,  high  frequency  decoupling  factors  in  the  30  to  50  Hz  frequency  band  are  larger  for 
overdecoupled  explosions  than  for  fully  decoupled  explosions.  Decoupling  factors  of  less 
than  one  are  predicted  for  frequencies  higher  than  about  100  Hz  for  explosions  in  tuff. 
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1.  INTRODUCTION 


It  is  well  known  that  detonating  an  explosion  in  a  cavity  decouples  the 
explosion  source  from  the  surrounding  material,  leading  to  sharply  reduced 
seismic  signals  compared  with  a  tamped  explosion  of  the  same  yield.  The 
possibility  of  such  explosions  being  used  as  a  means  of  avoiding  detection  has 
been  a  matter  of  serious  concern  for  treaty  verification.  The  theoretical 
decoupling  factor  is  a  strong  function  of  frequency,  and  is  largest  at  low 
frequencies.  At  high  frequencies  both  the  decoupling  factor  and  the  amplitude 
of  the  signal  are  reduced.  This  has  led  Evernden,  at  at.  (1986)  and  others  to 
suggest  that  it  may  be  possible  to  detect  and  identify  decoupled  explosions 
using  high  frequency  ( >  20  Hz)  seismic  data. 


The  purpose  of  this  study  is  to  gain  an  improved  understanding  of  the 
seismic  source  function  of  decoupled,  partially  coupled,  and  overdecoupled 
explosions  at  low  and  high  frequencies,  and  to  evaluate  the  implications  of  these 
predictions  for  detection  and  discrimination.  In  an  earlier  report  (Murphy,  et  at., 
1988),  we  performed  finite  difference  calculations  in  tuff  and  salt,  and  compared 
the  results  with  data  for  the  Mill  Yard  and  Sterling  experiments.  The  principal 
conclusions  of  that  report  based  on  this  preliminary  set  of  calculations  were  that: 

1.  A  step  pressure  response  is  an  inadequate  model  for  a 
decoupled  explosion.  The  air  shock  produces  enhanced 
seismic  coupling  at  high  frequencies  relative  to  a  step 
pressure  source. 

2.  The  Latter  decoupling  criterion  is  approximately  correct  for 
salt,  however  nonlinearities  due  to  pore  crushing  in  tuff 
may  occur  at  yields  substantially  less  than  the  Latter 
criterion. 

3.  Because  of  pore  crushing,  decoupling  effectiveness  in  tuff 
may  be  greater  than  decoupling  effectiveness  for  an 
equivalent  linear  material.  In  salt,  however,  where  the 
primary  nonlinear  effect  is  yielding,  the  nonlinear  effects 
lead  to  increased  coupling  effectiveness. 
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4.  Examination  of  tree-field  data  from  Sterling  and  Mill  Yard 
indicates  that  high  frequencies  couple  effectively  into  the 
seismic  regime. 

5.  The  observed  complexity  of  high  frequency  propagation 
effects  of  the  type  observed  at  ground  zero  above  the  Mill 
Yard  explosion  makes  it  unlikely  that  differences  in  high 
frequency  source  characteristics  can  be  reliably  recovered 
from  recorded  seismic  data. 

In  this  report,  we  extend  the  earlier  simulations  to  encompass  a  wider 
range  of  yields,  and  examine  the  seismic  source  functions  for  overdecoupled 
and  underdecoupled  explosions.  In  Sections  3  and  4  of  this  report,  we  present 
the  results  from  a  series  of  one  dimensional  finite  difference  simulations  covering 
a  range  of  yields  in  unsaturated  tuff  and  salt  cavities.  Each  set  of  calculations 
was  repeated  four  times:  first  with  a  full  nonlinear  model  of  the  cavity  airshock 
and  the  emplacement  material;  second  with  a  nonlinear  model  of  the  airshock. 
but  with  a  linear  material  external  to  the  cavity;  third  with  a  full  nonlinear  model  of 
the  material  outside  the  cavity,  but  with  the  explosion  applied  as  a  step  in 
pressure  on  the  cavity  wall;  and  fourth  with  a  step  pressure  applied  to  a  linear 
material.  The  linear  material  in  each  case  has  the  same  properties  as  the 
nonlinear  material,  except  that  it  does  not  yield  or  undergo  pore  crushing.  This 
set  of  calculations  allows  us  to  compare  the  full  nonlinear  results  with  results 
obtained  from  simpler  calculations.  In  addition,  we  have  performed  finite 
difference  simulations  of  tamped  explosions  using  the  full  nonlinear  models  for 
salt  and  tuff.  We  compare  these  results  with  the  decoupled  simulations  to 
determine  the  frequency  dependent  effects  of  decoupling. 
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2.  REVIEW  OF  CAVITY  DECOUPLED  EXPLOSION  THEORY 


In  this  section,  we  review  some  fundamental  theoretical  relations  for 
explosions  in  a  cavity.  These  relations  are  used  in  the  following  sections  of  this 
report  for  comparison  with  the  detailed  results  of  finite  difference  calculations. 
The  reduced  velocity  potential  for  a  pressure  pulse  P(t).  with  derivative  P(t),  and 
corresponding  Fourier  transforms  P(w)  and  P(w)  is  given  by: 


*(«)  =  P(w) 


q3  2 

R  wQ 
4/i 


iu/u 


0 


uo  "  iwow 


(X  +  2/0/4/!)],/ 


(2.1) 


where  =  a/R,  R  is  the  cavity  radius,  a  is  the  compressional  velocity  of  the 
external  medium,  and  X  and  ft  are  the  Lame  constants  of  the  external  medium. 
For  a  step  in  pressure  of  magnitude  PQ  applied  at  time  t  =  0,  P(w)  =  PQ.  and  for 
this  or  any  pressure  pulse  P(t)  with  static  value  PQ,  in  the  low  frequency  limit  i{u) 
becomes: 

=  P0  r3/4^  (2'2> 

For  an  explosion  in  an  air  filled  cavity,  the  static  value  of  the  pressure  is 
related  by  conservation  of  energy  to  the  explosion  yield  W  by: 


(l-U  w. 

V 


(2.3) 


3 

where  V  =  4/3tR  is  the  cavity  volume,  and  7  is  the  adiabatic  expansion 
constant  which  is  approximately  1.2  for  air.  From  Equations  2.2  and  2.3, 


*_  = 


16  *71 


(7-1)  W 


(2.4) 
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Note  that  for  a  decoupled  explosion,  ^  depends  quite  strongly  on  the 
shear  modulus  of  the  external  medium.  For  the  material  models  used  in  the 

q  g 

following  sections,  the  shear  moduli  are  3.5  x  10  Pascals,  and  14.2  x  10 

1 2 

Pascals,  for  tuff  and  salt,  respectively.  Using  a  conversion  factor  of  4.1  x  10 
Joules/Kiloton  for  a  nuclear  explosion,  we  find  for  a  decoupled  explosion  in  tuff 

*  ~  14.0  W  (2.5) 

00  ~ 

and  for  a  decoupled  explosion  in  salt 

t  ~  3.45  W  (2.6) 

3 

where  W  is  the  yield  in  kilotons  and  is  in  m  Variations  in  7  also  can  cause 
significant  variations  in 

The  relative  body  wave  amplitude  for  explosion  sources  in  dissimilar 
materials  with  the  same  propagation  path  is  proportional  to  fpa  #(w)  (Stevens 
and  Day,  1985),  where  p  and  a  are  the  density  and  compressional  velocity  of  the 
external  medium,  so  from  Equation  2.4  Ab  is  proportional  to  Vpixip  W  in  the  low 
frequency  limit.  Neglecting  higher  frequency  spectral  differences,  therefore.  A^ 
will  be  relatively  larger  for  decoupled  explosions  in  lower  velocity  media.  This 
difference  is  approximately  a  factor  of  three  for  the  salt  and  tuff  models 
considered  here. 

The  surface  wave  amplitude  Ag  of  a  decoupled  explosion  is  proportional 
to  jii(u),  so  Ag  is  independent  of  the  material  external  to  a  decoupled  explosion 
in  the  low  frequency  limit. 

In  the  high  frequency  limit,  Equation  2.1  reduces  to 

i  ’f  ■  ur 

Re  ^ 

¥  u  ~  -P(tf)  (2  7) 
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For  a  step  function  response,  P(«)  is  equal  to  the  static  pressure  PQ.  As  we 

discuss  m  the  following  two  sections,  however,  P(w)  may  be  quite  complex  at 

high  frequencies  because  of  cavity  reverberations  and  the  effect  of  the  shock 

wave  on  the  cavity  wall.  Also,  certain  nonlinear  effects  such  as  pore  crushing 

-2 

cause  the  spectrum  to  decay  more  rapidly  than  w  for  partially  coupled 
explosions. 

In  the  high  frequency  limit,  A^  becomes 


Ab  -  P(W) 


-  alp  R  u 


-2 


(2.8) 


so  the  dependence  on  material  properties  is  in  the  factor  fop.  High  frequency 
body  waves  are  therefore  less  sensitive  to  material  properties  than  lower 
frequency  body  waves.  For  the  material  models  considered  here,  this  factor  is  30 
percent  larger  for  salt  than  for  tuff. 

Commonly  used  criteria  for  full  decoupling  are  given  by  the  relation 


i  kpgh  ( 2.9) 


where  k  is  between  0,5  (Latter  criterion,  Latter,  ef  a/.,  1961),  and  1.0  (Patterson 
criterion,  Patterson.  1964),  and  ^gh  is  the  overburden  pressure  at  depth  h. 

In  the  following  two  sections,  we  consider  explosions  in  an  1 1  meter 
radius  tuff  cavity  with  62  bar  overburden  pressure,  and  a  17  meter  salt  cavity 
with  175  bar  overburden  pressure.  For  these  two  cases,  the  decoupling  criteria 
range  from  21  to  42  tons  for  tuff,  and  from  220  to  440  tons  for  salt. 
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3.  FINITE  DIFFERENCE  SIMULATIONS  OF  EXPLOSIONS 
IN  TUFF  CAVITIES 


We  performed  the  set  of  cavity  decoupled  simulations  described  in  the 
introduction  over  a  range  of  yields  from  5  to  160  tons  for  an  1 1  meter  radius 
cavity  in  unsaturated  tuff.  The  tuff  model  used  here  was  described  in  our  earlier 
report  (Murphy,  eta/.,  1988),  and  is  appropriate  for  Ranier  Mesa  tuff.  This  model 
has  an  initial  porosity  of  1 ,8  percent,  and  P  and  S  wave  speeds  of  2558  m'sec 
and  1372  m/sec  respectively. 


A  simulation  of  a  tamped  explosion  in  this  material  was  run  for 
comparison.  In  the  low  frequency  limit,  the  tamped  simulation  leads  to  the  result 


¥ 


tamped  ~ 


325  W 


(3.1) 


3 

where  W  is  the  yield  in  kilotons  and  ^tampecl  is  in  m  .  From  Equation  2.5, 
therefore,  the  theoretical  low  frequency  decoupling  factor  for  this  material, 
assuming  7  of  1.2  and  linear  behavior  of  the  medium ,  is  given  by: 


D  = 


ll|  /  lif  r* 

tamped  » ~ 


23 


(3.2) 


As  mentioned  in  the  previous  section,  the  decoupling  factor  is  a  material 

dependent  quantity,  and  it  depends  quite  strongly  on  both  the  shear  modulus  of 

the  material  and  the  tamped  coupling  of  the  material.  For  the  tamped  explosion, 

1  /3 

the  final  cavity  radius  is  14.3  m/KT  ,  the  radius  of  plastic  yielding  is 
104  m/KT1^3,  and  the  pore  crush  radius  is  157  m/KT1'3.  Note  that  in  this 
material  the  transition  to  elastic  behavior  occurs  at  the  limit  of  pore  crushing 
rather  than  at  the  limit  of  plastic  yielding. 


In  Tables  3.1  through  3.4,  we  summarize  the  results  of  these 
calculations.  The  columns  in  each  table  give  the  yield  of  the  simulation.  4'o>,  the 
peak  value  of  the  time  domain  reduced  displacement  potential  *DeaK,  the 
decoupling  ratio  D  (*tampec/*„)>  the  radius  to  limit  of  plastic  yielding  R  and  the 
radius  to  the  limit  of  pore  crushing  R 
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TABLE  3.3 


EXPLOSION  IN  11  METER  TUFF  CAVITY  -  STEP  PRESSURE 


Yield 

kt 


Vpegk 

mi 


0.040 

0.477 

0.540 

27.3 

11.00 

11.00 

0.060 

0.689 

0.781 

28.3 

11.00 

11.00 

0.080 

0.940 

1.067 

27.7 

11.00 

11.00 

0.160 

2.009 

2.275 

25.9 

11.00 

11.00 

TABLE  3.4 

EXPLOSION  IN  11  METER  TUFF  CAVITY  - 
STEP  PRESSURE:  LINEAR 


Yield 

kt 


11.30 

11.00 

11.30 

11.00 

11.30 

11.00 

11.30 

11.00 

11.30 

12.32 

11.30 

16.57 

16.93 

29.87 

SSS-TR-BS-lGoS 


In  Figure  3.1,  we  show  the  low  frequency  decoupling  factor  for  the  air- 
filled  linear  and  nonlinear  cases  listed  in  Tables  3.1  and  3.2.  The  Latter  criterion 
for  full  decoupling  in  this  material  is  21  tons.  As  can  be  seen  from  the  tables,  the 
nonlinear  material  behavior  is  too  complex  to  be  described  by  a  single  number. 
For  a  model  with  a  true  air  shock  inside  the  cavity,  yielding  begins  with  an 
explosion  yield  of  between  20  and  40  tons,  however  pore  crushup  begins  with 
explosion  yields  as  small  as  5  tons.  The  effect  of  pore  crushup  is  to  remove 
energy  from  the  outgoing  seismic  wave,  so  it  has  the  effect  of  increasing 
decoupling.  As  a  result,  the  maximum  low  frequency  decoupling  factor  of  44  is 
achieved  with  a  yield  of  20  tons,  which  exceeds  the  decoupling  factor  of  23 
predicted  with  a  linear  material  model.  With  higher  explosion  yields,  the 
decoupling  effectiveness  decreases  as  yielding  becomes  the  dominant  effect.  At 
160  tons,  the  decoupling  factor  decreases  to  less  than  that  for  a  linear  model. 
For  the  linear  models,  the  low  frequency  decoupling  ratio  is  close  to  the  value  of 
23  predicted  from  simple  theory.  The  differences  from  this  value  result  from  the 
differences  in  the  equation  of  state  for  air  at  different  energy  levels. 

In  Figures  3.2  through  3.6,  we  show  the  reduced  velocity  potentials  for 
all  seven  yields  for  each  set  of  calculations.  Care  was  taken  in  the  calculations 
to  preserve  the  high  frequencies  by  using  fine  zoning  and  a  low  value  for  artificial 
viscosity.  The  results  are  accurate  to  approximately  500  hz.  Above  this 
frequency,  the  spectral  amplitudes  are  gradually  reduced  by  artificial  viscosity. 
All  of  the  cavity  calculations  containing  an  airshock  show  a  high  frequency  peak 
corresponding  to  cavity  reverberations.  For  the  smallest  yield  of  5  tons,  the  air 
shock  reverberations  dominate  the  spectrum.  The  dominant  frequency  of  this 
peak  increases  with  yield,  ranging  from  100  Hz  for  the  5  ton  calculation  to  600  hz 
for  the  160  ton  calculation.  This  result  indicates  that  the  high  frequency 
enhancement  due  to  cavity  reverberations  is  more  likely  to  be  observed  from  an 
overdecoupled  explosion  than  from  a  fully  decoupled  or  partially  decoupled 
explosion.  Note  that  these  frequencies  will  cube  root  scale  to  other  yields  and 
geometries  provided  that  the  ratio  of  the  yield  to  cavity  volume  is  kept  constant. 
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Figure  3.1.  Low  frequency  decoupling  factor  plotted  as  a  function  of  yield 
for  the  air-filled  linear  and  nonlinear  tuff  simulations.  The 
decoupling  factor  is  defined  to  be  the  ratio  of  the  tamped 
reduced  velocity  potential  to  the  decoupled  reduced  velocity 
potential  with  the  same  material  properties  used  for  both 
calculations.  The  vertical  line  indicates  the  yield  corresponding 
to  the  Latter  decoupling  criterion  for  this  material  and  geometry. 
The  peak  in  the  nonlinear  decoupling  factor  is  due  to  the  effect 
of  pore  crushing  which  reduces  the  low  frequency  source 
function. 
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Figure  3.4. 
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Reduced  velocity  potential  for  seven  different  yields  for  an  air 
filled  cavity  inside  a  medium  with  a  linear  tuff  material  model. 
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Figure  3.5.  Reduced  velocity  potential  for  seven  different  yields  for  a  step 
pressure  on  a  cavity  inside  a  medium  with  a  nonlinear  tuff 
material  model. 


Figure  3.6. 


Reduced  velocity  potential  for  seven  different  yields  for  a  step 
pressure  on  a  cavity  inside  a  medium  with  a  linear  tuff  material 
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In  Figures  3.7  through  3.10,  we  show  the  decoupling  factors  for  each  of 
the  simulations,  plotted  as  a  function  of  frequency.  These  decoupling  factors 
were  determined  by  cube-root  scaling  the  tamped  explosion  spectrum  to  the 
same  yield  as  the  cavity  simulation,  and  then  taking  the  ratio  of  tamped  to  cavity 
spectra.  At  low  frequencies,  there  are  significant  differences  between  the  linear 
and  nonlinear  simulations.  At  high  frequencies,  however,  the  cavity 
reverberations  dominate  the  spectrum.  The  only  significant  difference  between 
the  linear  and  nonlinear  calculations  at  high  frequency  is  that  the  nonlinear 
spectra  are  attenuated  by  pore  crushing.  Because  the  cavity  reverberations 
generate  a  substantial  amount  of  high  frequency  energy,  the  high  frequency 
decoupling  factors  are  significantly  lower  for  the  air-filled  cavity  simulations  than 
for  the  step  pressure  calculations.  Because  the  nonlinear  effects  cause 
attenuation  of  high  frequencies  in  the  tamped  case,  however,  the  decoupling 
factors  drop  off  rapidly  at  high  frequency,  and  are  reduced  to  less  than  one  at 
sufficiently  high  frequency  in  all  cases. 


16 


Decoupling  Factor 


SSS-TR-89- 10755 


160  Ton 
80  Ton 
60  Ton 
40  Ton 
20  Ton 
10  Ton 
5  Ton 


I)  f 


Frequency  (Hz) 


Figure  3.7.  Frequency  dependent  decoupling  factor  for  seven  different 
yields  for  an  air  filled  cavity  inside  a  medium  with  a  nonlinear  tuff 
material  model. 
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Figure  3.8. 


Frequency  dependent  decoupling  factor  for  seven  different 
yields  for  an  air  filled  cavity  inside  a  medium  with  a  linear  tuff 
material  model. 
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Figure  3.9. 


Frequency  dependent  decoupling  factor  for  seven  different 
yields  for  a  step  pressure  on  a  cavity  inside  a  medium  with  a 
nonlinear  tuff  material  model. 
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4.  FINITE  DIFFERENCE  SIMULATIONS  OF 
EXPLOSIONS  IN  SALT  CAVITIES 

We  performed  the  set  of  cavity  decoupled  simulations  described  in  the 
introduction  over  a  range  of  yields  from  95  to  1520  tons  for  a  17  meter  radius 
cavity  m  salt.  The  constitutive  model  used  in  these  calculations  was  based  on 
.ne  results  of  laboratory  material  property  tests  on  polycrystalline  salt  (Heard, 
et  a!..  1975).  We  also  made  some  test  calculations  using  the  earlier 
hardening/softening  model  of  Rimer  and  Cherry  (1982).  This  earlier  model 
produced  an  excellent  comparison  with  the  nearfield  Salmon  data,  however  it 
was  found  to  be  too  weak  to  maintain  a  17  meter  cavity  at  this  depth  without 
undergoing  very  substantial  nonlinear  deformation.  A  test  calculation  using  a  17 
meter  cavity  in  this  material  showed  that  it  stabilizes  after  undergoing  a  surface 
displacement  (inward)  of  2  centimeters,  and  plastic  deformation  extending  to  70 
meters.  Since  the  Sterling  cavity  was  generated  by  the  Salmon  explosion, 
however,  it  is  not  clear  what  initial  conditions  to  use  for  studying  partially  coupled 
explosions  in  this  material.  A  cavity  constructed  in  this  material  would  develop  a 
hardened  core  extending  to  a  considerable  distance  away  from  the  cavity.  The 
calculations  in  this  section  have  been  performed  using  the  polycrystalline  salt 
model  described  above  which  is  strong  enough  to  support  a  stable  cavity. 
Some  calculations  have  also  been  performed  using  the  much  weaker 
Rimer/Cherry  model  for  the  air-filled,  partially  coupled  cases,  and  the  results  are 
shown  for  comparison  in  some  of  the  figures  in  this  section. 

In  Figure  4.1,  we  show  a  comparison  of  particle  velocities  at  four 
locations  from  the  simulation  at  0.38  KT,  with  the  observed  particle  velocities 
from  the  Sterling  explosion.  The  data  cannot  be  fit  exactly  with  a  spherically 
symmetric  model,  since  measured  particle  velocities  at  the  same  radial  distance 
at  angles  obliquely  above  and  below  the  explosion  differ  in  amplitude  by  about  a 
factor  of  two.  Nevertheless,  the  model  appears  to  provide  a  reasonably  good 
representation  of  the  Sterling  data. 
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Figure  4.1 .  Comparison  of  data  from  the  Sterling  explosion  with  calculation 
of  a  0.38  KT  explosion  in  an  air  filled  17  meter  cavity  using  the 
nonlinear  material  model  described  in  the  text.  The  data  points 
were  advanced  by  6  msec.  14B-27  is  at  shot  depth,  166  meters 
from  the  explosion.  14C-22  is  above  shot  depth,  225  meters 
from  the  explosion.  14C-32  is  below  shot  depth,  225  meters 
from  the  explosion.  6-27  is  at  shot  depth,  318  meters  from  the 
explosion. 
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A  simulation  of  a  tamped  explosion  in  this  material  was  run  for 
comparison.  In  the  low  frequency  limit,  the  tamped  simulation  gives  the  result 


I' 

tamoea 


=  320  W 


(4.1) 


where  W  is  the  yield  in  kilotons.  From  Equation  2.6,  therefore,  the  theoretical  low 
frequency  decoupling  factor  for  this  material  is  given  by: 


D  = 


iii 

tamped 


~  93 


(4.2) 


In  Tables  4.1  through  4.4,  we  summarize  the  results  of  these 

calculations.  The  columns  in  each  table  give  the  yield  of  the  simulation,  4^,  the 

peak  value  of  the  time  domain  reduced  displacement  potential  4>peak,  the 

decoupling  ratio  D  (*tampecl/*J.  and  the  radius  to  limit  of  plastic  yielding  R 

There  is  no  pore  crush  radius  in  this  case,  because  the  salt  is  assumed  to  be 

1 

nonporous.  For  this  tamped  explosion,  the  cavity  radius  is  12.0  m/KT  ,  and 
the  elastic  radius  is  112  m/KT1^3 


In  Figure  4.2,  we  show  the  low  frequency  decoupling  ratio  plotted  as  a 
function  of  yield  for  the  air-filled  linear  and  nonlinear  cases.  We  also  show 
results  calculated  using  the  Rimer/Cherry  model  for  comparison.  The  Latter 
criterion  for  full  decoupling  in  this  material  is  220  tons.  As  can  be  seen  from  the 
tables  and  Figure  4.2,  yielding  begins  at  slightly  less  than  190  tons  for  the  full 
nonlinear  case.  With  the  step  function  source  model,  there  is  no  plastic  yielding 
for  explosion  yields  less  than  about  a  kiloton  for  this  material.  The  decoupling 
factor,  however,  remains  about  the  same  over  a  large  range  of  explosion  yields, 
and  decreases  by  only  a  factor  of  2  at  1520  tons,  even  though  the  elastic  radius 
has  increased  to  33  meters.  This  slow  decrease  in  decoupling  is  also  true  for 
the  Rimer/Cherry  model.  Even  though  this  model  is  initially  very  weak,  and  is 
not  linear  at  any  yield  used  in  this  set  of  simulations,  the  decoupling  factor  has 
decreased  by  only  about  a  factor  of  four  at  1520  tons  yield.  It  can  be  seen  from 
this  figure  that  the  main  features  of  the  decoupling  calculations  are  not 
particularly  sensitive  to  the  material  model  employed  in  that  all  of  the 
calculations  indicate  that  the  low  frequency  decoupling  effectiveness  decreases 
quite  slowly  with  increasing  yield  above  the  Latter  threshold. 
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TABLE  4.1 

EXPLOSION  IN  17  METER  AIR-FILLED  SALT  CAVITY 


Yield 

kt 

m3 

Vpeqk 

mi 

D 

R, 

m 

0.095 

0.297 

0.413 

101.7 

17.00 

0.138 

0.453 

0.613 

96.9 

17.00 

0.190 

0.732 

0.904 

82.6 

17.94 

0.276 

1.151 

1.397 

76.3 

18.53 

0.380 

1.683 

2.003 

71.8 

19.49 

0.570 

2.909 

3.339 

62.3 

20.82 

0.760 

4.135 

4.723 

58.5 

22.01 

1.520 

10.375 

11.545 

46.6 

33.27 

TABLE  4.2 

EXPLOSION  IN  1 7  METER  AIR-FILLED  SALT  CAVITY:  LINEAR 


Yield 

kt 

Voo 

m3 

Vpeqk 
m 3 

D 

Ry 

m 

0.297 

0.413 

101.7 

— 

0.453 

0.613 

96.9 

0.190 

0.721 

0.860 

83.8 

17.00 

0.276 

0.996 

1.207 

88.1 

17.00 

0.380 

1.309 

1.560 

92.3 

17.00 

0.570 

1.962 

2.259 

92.4 

17.00 

0.760 

2.597 

2.989 

93.1 

17.00 

1.520 

5.669 

6.367 

85.3 

17.00 
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TABLE  4.3 

EXPLOSION  IN  17  METER  SALT  CAVITY  -  STEP  PRESSURE 


Yield 

kt 

m3 

Vpeak 

mi 

D 

Ry 

m 

0.095 

0.294 

0.340 

102.8 

17.00 

0.138 

0.457 

106.3 

17.00 

0.606 

111.5 

17.00 

0.774 

0.889 

113.4 

17.00 

1.149 

mmm 

105.2 

17.00 

1.753 

mEm 

103.4 

17.00 

2.332 

2.615 

103.7 

17.00 

6.321 

6.859 

76.5 

24.61 

TABLE  4.4 

EXPLOSION  IN  17  METER  SALT  CAVITY  - 
STEP  PRESSURE:  LINEAR 


Yield 

kt 

v« 

_  3 

m 

Vpeqk 

mi 

D 

Ry 

m 

0.294 

0.340 

102.8 

17.00 

0.138 

0.413 

0.457 

106.3 

17.00 

0.542 

0.606 

111.5 

17.00 

0.276 

0.774 

0.889 

113.4 

17.00 

*1’"“-’  \wis 

1.149 

1.280 

105.2 

17.00 

1 

1.753 

1.970 

103.4 

17.00 

2.332 

2.615 

103.7 

17.00 

5.241 

5.872 

92.2 

17.00 
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Figure  4.2.  Low  frequency  decoupling  factor  plotted  as  a  function  of  yield 
for  the  air-filled  linear  and  nonlinear  salt  simulations.  Also  shown 
are  decoupling  factors  calculated  using  the  weaker 
Rimer/Cherry  salt  model.  The  vertical  line  indicates  the  yield 
corresponding  to  the  Latter  decoupling  criterion  for  this  material 
and  geometry.  The  decoupling  factors  decrease  by  only  a 
factor  of  between  two  and  four  at  the  highest  yield  calculated. 
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In  Figure  4.3  we  show  a  comparison  of  the  normalized  low  frequency 
decoupling  factors  obta med  from  our  calculations  and  those  of  Patterson  (1964) 
with  the  observed  Decoupling  values  from  the  COWBOY  HE  decoupling  tests 
which  covered  the  widest  available  range  of  ratios  of  yield  to  cavity  volume.  In 
order  to  compare  the  dependence  of  decoupling  effectiveness  on  cavity  size,  the 
three  theoretical  curves  in  Figure  4.3  have  been  normalized  to  the  observed 
COWBOY  full  decoupling  value  of  70.  It  can  be  seen  that  all  three  sets  of 
calculations  are  in  qualitative  agreement  with  the  limited  data  and  indicate  that 
the  decoupling  effectiveness  decreases  fairly  slowly  with  increasing  yield  over 
the  range  investigated  in  this  report.  However,  in  detail,  the  somewhat  more 
rapid  decrease  predicted  using  the  weak  Rimer/Cherry  salt  model  does  seem  to 
provide  somewhat  better  agreement  with  the  two  COWBOY  HE  data  points 
associated  with  the  largest  overdrive  values. 

In  Figures  4.4  through  4.8,  we  show  the  reduced  velocity  potentials  for 
all  eight  yields  for  each  set  of  calculations.  As  with  the  tuff  models,  all  of  the 
calculations  containing  an  airshock  show  a  high  frequency  peak  corresponding 
to  cavity  reverberations.  Again,  the  frequency  of  the  peak  increases  with  yield, 
ranging  from  just  over  100  Hz  for  the  95  ton  simulation,  to  about  600  Hz  for  the 
1520  ton  simulation.  In  all  cases,  there  is  little  difference  between  the  linear  and 
nonlinear  simulations  at  high  frequencies. 

In  Figures  4.9  through  4.12,  we  show  the  decoupling  ratio  for  each  of  the 
simulations.  The  results  are  very  similar  to  the  tuff  simulations.  The  air-filled 
cavity  simulations  have  lower  decoupling  ratios  than  the  step  pressure 
simulations  at  high  frequencies,  particularly  for  overdecoupled  explosions. 
Decoupling  ratios  drop  off  at  high  frequencies,  however  the  dropoff  occurs  at 
higher  frequencies  for  lower  yield  explosions  than  for  higher  yield  explosions. 
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Figure  4.3.  Low  frequency  decoupling  factors  for  the  air-filled  nonlinear  salt 
model,  together  with  decoupling  factors  calculated  with  the 
Rimer/Cherry  model,  the  finite  difference  calculations  of 
Patterson  (1964),  and  data  points  from  the  Cowboy  experiment. 
These  are  plotted  against  the  dimensionless  quantity  r/rd  where 
rd  is  the  radius  for  “full"  decoupling. 
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Reduced  velocity  potential  for  a  tamped  explosion  in  salt,  scaled 
to  eight  different  yields. 
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Figure  4.6.  Reduced  velocity  potential  for  eight  different  yields  for  an  air 
filled  cavity  inside  a  medium  with  a  linear  salt  material  model. 
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Figure  4.7. 


Reduced  velocity  potential  for  eight  different  yields  for  a  step 
pressure  on  a  cavity  inside  a  medium  with  a  nonlinear  salt 
material  model. 
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Figure  4.8. 


Reduced  velocity  potential  for  eight  different  yields  for  a  step 
pressure  on  a  cavity  inside  a  medium  with  a  linear  salt  material 
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Figure  4.9. 


Frequency  dependent  decoupling  factor  for  eight  different  yields 
for  an  air  filled  cavity  inside  a  medium  with  a  nonlinear  salt 
material  model. 
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Figure  4.10.  Frequency  dependent  decoupling  factor  for  eight  different  yields 
for  an  air  filled  cavity  inside  a  medium  with  a  linear  salt  material 
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Figure  4.1 1 .  Frequency  dependent  decoupling  factor  for  eight  different  yields 
for  a  step  pressure  on  a  cavity  inside  a  medium  with  a  nonlinear 
salt  material  model. 
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Figure  4.12.  Frequency  dependent  decoupling  factor  for  eight  different  yields 
for  a  step  pressure  on  a  cavity  inside  a  medium  with  a  linear  salt 
material  model. 
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5.  DISCUSSION 

An  important  result  of  this  study  is  that  we  predict  that  a  substantial 
amount  of  decoupling  can  be  obtained  using  cavities  that  are  significantly 
smaller  than  would  be  required  for  full  decoupling.  This  result  is  consistent  with 
earlier  studies.  For  example,  the  OTA  report  (Office  of  Technology  Assessment, 
1988,  Figure  6-4),  shows  that  a  10  meter  radius  cavity  is  sufficient  to  decouple  a 
one  kiloton  explosion  to  within  about  a  factor  of  two  of  full  decoupling,  even 
though  full  decoupling  would  require  a  25  meter  radius  cavity. 

In  this  section  we  discuss  three  related  topics:  the  imp'ications  of  this 
result  for  detection  of  decoupled  explosions;  the  possibility  of  discriminating 
decoupled  explosions  from  other  seismic  sources  because  of  spectral 
differences;  and  the  stability  of  cavities  used  for  detonating  partially  coupled 
explosions. 

Rewriting  the  tables  in  the  previous  sections  as  a  function  of  the  ratio  of 
yield  W  to  the  yield  required  for  full  decoupling  Wf,  we  find  the  following  low 
frequency  decoupling  factors  (D)  listed  in  Tables  5.1  and  5.2. 

The  significance  of  these  results  for  evasion  assessment  can  be 
envisioned  more  easily  by  considering  mb  values  which  would  be  expected  from 
explosions  of  different  yields  in  cavities  which  span  the  range  of  sizes  generally 
thought  to  be  feasible  (see,  for  example  the  1988  OTA  report).  For  purposes  of 
illustration,  we  wi<l  assume  that  coupling  in  salt  is  comparable  to  that  in 
competent  rock  media  typified  by  that  below  the  water  table  at  NTS.  An 
appropriate  mb/yield  curve  for  such  media  has  been  reported  by  Murphy  (1981) 
to  be 


mb  =  3.94  +  0.81  log  W  .  (5.1) 

Similarly,  it  has  been  observed  at  NTS  that  for  explosions  in  dry  tuff,  the 
observed  mb  value  for  a  given  yield  is  lower  than  predicted  by  equation  5.1  by 
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TABLE  5.1 

DECOUPLING  FACTORS  IN  DRY  TUFF 


wtwf 

D 

log  D 

0.95 

44 

1.64 

1.91 

36 

1.55 

2.86 

32 

1.50 

3.81 

29 

1.45 

7.62 

19 

1.27 

TABLE  5.2 

DECOUPLING  FACTORS  IN  SALT 


W/Wf 

D 

0.86 

83 

1.92 

1.25 

76 

1.88 

1.73 

72 

1.86 

2.59 

62 

1.79 

3.45 

59 

1.77 

6.91 

47 

1.67 
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0.5  to  1.0  magnitude  units  (Murphy,  1981).  Taking  0.6  as  a  representative  value, 
we  obtain  a  corresponding  NTS  m^/yield  curve  for  dry  tuff  of  the  form 

mb  =  3.34  -  0.81  log  W  .  (5.2) 

For  explosions  in  the  Soviet  Union,  it  is  generally  believed  that  the  mb 
value  for  a  given  yield  can  be  expected  to  be  between  0.3  and  0.4  magnitude 
units  greater  than  those  given  by  Equations  5.1  and  5.2  (OTA.  1988).  Adopting 
0.3  as  our  nominal  bias  value,  we  then  have  for  tamped  explosions  in  Soviet  salt 

mb  =  4.24  -  0.81  log  W  (5.3) 

and  for  tamped  explosions  in  dry  tuff 

mb  =  3.64  -  0.81  log  W  .  (5.4) 

Making  the  reasonable  assumption  that  the  dependence  of  decoupling 
effectiveness  on  yield  is  independent  of  the  absolute  value  of  (i.  e.  simple 
yield  scaling  is  applicable),  it  is  a  straightforward  matter  to  estimate  the  mb/yield 
relations  corresponding  to  explosions  in  different  sized  cavities.  Thus,  if  mb  (T) 
denotes  the  mb  values  corresponding  to  a  fully  tamped  explosion  of  a  given 
yield  and  mb  (D)  denotes  the  mb  value  corresponding  to  that  same  yield 
explosion  in  a  particular  sized  cavity,  then 

mb  (D)  =  mb  (T)  -  log  D  (5.5) 

where  D  is  the  decoupling  factor  appropriate  for  that  value  of  W./W^. 

In  Figures  5  1  and  5.2,  we  show  the  results  of  evaluating  Equation  5.5  for 
cavities  large  enough  to  fully  decouple  1  kiloton,  5  kiloton,  and  10  kiloton 
explosions  according  to  the  Latter  criterion.  The  figures  show  that  if  we  adhere 
to  the  mb  3.5  identification  threshold  quoted  in  the  1988  OTA  report,  decoupled 
explosions  with  yields  on  the  order  of  20  to  30  kilotons  could  evade 
identification,  even  if  it  were  only  possible  to  construct  a  cavity  large  enough  to 
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Predicted  body  wave  magnitude  for  explosions  in  dry  tuff.  The 
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Predicted  body  wave  magnitude  for  explosions  in  salt.  The  four 
curves  indicate  the  mb  for  a  tamped  explosion,  and  m^  for 
explosions  in  cavities  sufficient  to  decouple  i  KT,  5  KT,  and 
10  KT  explosions. 
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fully  decouple  5  KT.  In  fact,  the  results  shown  in  these  figures  suggest  that  it 
should  be  possible  to  detonate  explosions  with  yields  on  the  order  of  10  KT  in  a 
cavity  which  is  only  large  enough  to  fully  decouple  1  KT  and  still  remain  well 
below  the  projected  3.5  identification  threshold. 

Moving  to  the  issue  of  discriminating  decoupled  explosions,  we  note  that 
the  spectra  of  decoupled  explosions  are  distinctly  different  from  tamped 
explosions  because  of  the  high  frequency  cavity  oscillation  peak.  If  this  could  be 
detected  it  might  provide  a  way  to  discriminate  a  decoupled  explosion  from  a 
small  tamped  explosion.  Unfortunately,  the  peak  occurs  at  very  high  frequency. 
In  the  examples  shown  in  this  report,  the  peaks  all  occur  at  frequencies  greater 
than  100  Hz.  The  peak  frequency  will  cube-root  scale  to  lower  frequencies  for 
larger  yields  and  cavity  sizes.  For  a  hypothetical  fully  decoupled  5  kiloton 
explosion  in  a  48  meter  salt  cavity,  the  spectral  peak  would  occur  at  about  75  hz. 
This  peak  would  be  expected  to  occur  at  somewhat  lower  frequencies  for 
overdecoupled  explosions  and  at  somewhat  higher  frequencies  for  partially 
decoupled  explosions.  However,  it  appears  that  in  most  cases  of  practical 
concern,  it  would  occur  in  the  high  frequency  band  in  which  it  has  been  found  to 
be  very  difficult  to  separate  source  and  propagation  path  effects  (Murphy,  et  al., 
1988).  It  therefore  seems  unlikely  that  such  spectral  peaks  could  be  reliably 
detected  and  identified  at  regional  distances. 

Finally,  we  want  to  mention  the  difficult  issue  of  cavity  stability.  This  is 
important  for  two  reasons.  First,  an  evader  would  have  to  be  concerned  about 
detection  if  a  crater  formed  or  radioactive  products  leaked  from  the  explosion. 
Second,  given  the  high  cost  of  constructing  such  cavities,  it  would  be 
advantageous  if  the  cavities  were  reusable.  There  are  a  few  observations  that 
shed  some  light  on  this  question.  It  is  clear  that  cavity  stability  is  strongly 
material  dependent.  Cavities  resulting  from  tamped  explosions  in  salt,  for 
example  (e.g.  Salmon  and  Gnome),  appear  to  be  stable  in  general.  On  the  other 
hand,  cavities  resulting  from  tamped  or  nearly  tamped  explosions  in  tuff  nearly 
always  collapse.  For  decoupled  explosions  in  tuff,  we  know  that  smaller  cavities 
are  more  likely  to  stand  than  larger  cavities.  Red  Hot,  which  was  substantially 
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overdriven  collapsed.  Diamond  Dust  and  Diamond  Mine,  however,  which  were 
overdriven  to  a  lesser  degree,  stayed  open  at  least  for  the  first  100  seconds  after 
the  explosion  (whether  or  not  they  remained  open  to  later  times  is  unknown). 
Fully  decoupled  explosions  in  tuff,  such  as  Mill  Yard  and  Mini  Jade,  appear  to  be 
stable. 


The  current  state  of  the  art  in  numerical  modeling  is  not  adequate  to 
definitively  answer  the  question  of  cavity  stability,  particularly  for  late  time 
collapses.  The  capability  does  exist,  however,  to  assess  the  danger  of 
hydrofracture,  which  is  the  most  serious  danger  during  and  soon  after  cavity 
expansion  and  rebound. 

The  possibility  of  hydrofracture  is  increased  when  an  event  is  partially 
decoupled  by  emplacement  in  a  large  underground  cavity.  Parametric  ground 
motion  calculations  (Rimer,  etal.,  1987)  have  shown  that,  for  a  tamped  or  nearly 
tamped  event,  cavity  expansion,  and  subsequent  dynamic  inward  motion  (cavity 
rebound)  which  occurs  within  the  first  second,  induce  compressive  residual 
hoop  stresses  in  the  rock  surrounding  the  cavity  which  are  much  greater  than 
the  pressure  of  the  radioactive  cavity  gases.  (These  stresses  have  been 
measured  on  both  HE  and  nuclear  events  in  tuff.)  Thus,  even  if  hydrofractures 
are  initiated  dynamically  from  the  cavity  during  its  outward  expansion,  the 
compressive  hoop  stresses  developed  during  cavity  rebound  prevent  these 
hydrofractures  from  propagating  more  than  about  one  cavity  radius  from  the  wall 
of  the  expanded  cavity.  For  somewhat  more  decoupled  events,  the  smaller 
initial  cavity  pressure  results  in  less  cavity  expansion  but  higher  later  time  cavity 
pressure  and  lower  residual  stresses  than  for  tamped  event,  thus  increasing  the 
hydrofracture  possibility.  Late  time  cavity  pressure  increases  as  decoupling  is 
increased  until  only  a  small  fraction  of  the  device  energy  goes  into  cavity 
expansion.  For  larger  initial  cavities,  residual  stresses  are  negligible,  and  initial 
cavity  pressure  decreases  with  initial  cavity  size.  For  a  fully  decoupled  event  the 
cavity  pressure  is  approximately  equal  to  the  in  situ  stress  so  that  the 
hydrofracture  threat  is  minimal. 
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S-CUBED  has  recently  developed  a  capability  to  simultaneously  model 
the  dynamic  ground  motion  and  hyorofracture  propagation.  This  capability 
(Rimer,  eta!..  1987).  (Nilson.  era/..  1987),  called  F-CUBED  (Fast,  Fluid,  Fracture) 
models  the  fluid  flow  in  the  hydrofracture  including  loss  mechanisms  such  as 
fluid  friction,  steam  condensation,  heat  conduction  and  seepage  into  the  walls  of 
the  fracture.  These  models  have  been  applied  to  DNA  nuclear  events  in  tuff 
cavities,  LANL  events  in  alluvium  cavities,  and  tamped  LLNL  events  in 
complicated  hard  rock  geologies  on  Pahute  Mesa.  Two  examples  of  interest 
here  are  0  5  and  2.0  kiloton  overburied  events  in  11.0m  radius  hemispheric 
cavities  in  a  Rainier  Mesa  saturated  zeolitic  tuff  geology.  These  events  were 
overdriven  by  approximately  factors  of  24  and  95,  respectively,  which  are 
considerably  larger  overdrives  than  were  considered  in  the  calculations  in  this 
report.  For  the  0.5  kiloton  case,  the  calculated  final  cavity  volume  was  about 
double  its  initial  volume  and  a  worst  case  single  hydrofracture  was  calculated  to 
propagate  210  m  before  stopping  (depth  of  burial  was  approximately  400  m). 
The  fracture  stopped  as  a  result  of  lowering  of  cavity  pressure  due  to  gas  flow 
into  the  crack.  If  a  discrete  number  of  fractures  out  of  the  cavity  were  presumed, 
the  fracture  lengths  would  be  significantly  less  due  to  greater  cavity  gas 
depletion.  For  the  2.0  KT  case,  the  calculated  final  cavity  volume  was  1 1  times 
initial  volume  and  the  fracture  propagated  281  m  before  ending  in  a  very  porous 
tuff  layer  above  the  working  point.  Although  compressive  residual  stresses  were 
calculated  in  this  case,  they  were  not  large  enough  to  prevent  fracture 
propagation.  These  results  are  material  dependent,  with  more  brittle,  stronger 
rock  expected  to  induce  greater  hydrofracturing. 

In  the  calculations  done  in  this  report,  the  final  displacement  of  the  cavity 
surface  was  0.3  percent  of  the  cavity  radius  for  salt  and  0.5  percent  of  the  cavity 
radius  in  tuff  for  the  largest  cases  of  overdrive  considered  (approximately  a 
factor  of  8),  and  the  final  cavity  volumes  were  increased  by  10  percent  and  15 
percent  respectively.  Since  the  danger  of  hydrofracture  would  be  considerably 
less  than  in  the  cases  considered  above,  it  would  appear  to  be  possible  to 
construct  cavities  that  could  be  overdriven  by  this  amount  in  which  the  danger  of 
hydrofracture,  at  least  to  the  extent  that  detection  would  be  a  danger,  could  be 
minimized. 
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6.  CONCLUSIONS 


We  have  performed  a  series  of  nonlinear  finite  difference  calculations  to 
model  explosions  in  air-filled  cavities  in  salt  and  tuff  media.  The  purpose  of  the 
study  was  to  obtain  an  improved  understanding  of  the  seismic  source  function  of 
decoupled,  overdecoupled,  and  partially  coupled  explosions  The  results  of  this 
set  of  simulations  indicate  that: 

1.  The  low  frequency  decoupling  level  is  strongly  material 
dependent,  and  is  a  function  of  the  shear  modulus  and 
tamped  coupling  of  the  medium  and  to  a  lesser  extent  a 
function  of  the  air  inside  the  cavity. 

2.  The  amount  of  decoupling  decreases  slowly  as  yield  is 
increased  above  the  yield  required  for  complete  decoupling. 

In  both  the  tuff  and  salt  models  used  in  this  study,  low 
frequency  decoupling  is  reduced  by  only  about  a  factor  of 
two  at  an  explosion  yield  eight  times  the  Latter  decoupling 
criterion.  In  fact,  assuming  that  simple  yield  scaling  is 
applicable,  these  calculational  results  suggest  that  it  should 
be  possible  to  decouple  explosions  with  yields  on  the  order 
of  20  to  30  kilotons  in  a  cavity  large  enough  to  fully  decouple 
only  5  kilotons  and  still  remain  below  the  mb  3.5  identification 
threshold  cited  in  the  recent  1988  OTA  report  on  seismic 
monitoring 

3.  Although  the  current  state-of-the-art  in  numerical  modeling  is 
not  adequate  to  definitively  address  the  question  of  cavity 
stability,  available  calculational  results  for  inouced  fracturing 
and  limited  test  data  do  not  seem  to  suggest  that  cavity 
stability  should  be  a  limiting  factor  for  the  range  of  cavity 
overdrive  factors  considered  in  this  report. 

4.  In  salt  cavities,  plastic  flow  begins  at  about  the  yeld  predicted 
by  the  Latter  criterion.  However,  as  mentioned  above,  the 
amount  of  decoupling  decreases  slowly  with  increasing 
explosion  yield. 
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5.  Explosions  in  cavities  in  unsaturated  tuff  decouple  more  than 
would  be  predicted  based  on  linear  theory.  This  happens 
because  the  dominant  nonlinear  mechanism  which  is 
operating  at  low  stress  levels  is  pore  crushing  which 
attenuates  the  outgoing  wave.  A  secondary  result  of  this  is 
that  in  tuff  the  decoupling  factor  decreases  to  less  than  one  at 
sufficiently  high  frequencies.  At  sufficiently  high  explosion 
yields,  about  a  factor  of  five  higher  than  the  Latter  criterion, 
plastic  flow  becomes  the  dominant  mechanism,  and 
decoupling  decreases  to  less  than  predicted  by  elastic 
theory. 

6.  The  decoupling  factor  decreases  at  high  frequencies, 
however  this  decrease  occurs  at  higher  frequency  for  lower 
yield  events.  Because  of  this,  high  frequency  decoupling 
factors  are  larger  for  overdecoupled  explosions  than  for 
"fully”  decoupled  explosions  at  the  same  frequency  with  the 
same  cavity  geometry. 

7.  At  high  frequencies,  there  is  little  difference  between  the 
source  functions  predicted  by  linear  theory  and  those 
predicted  by  a  nonlinear  external  material  model,  except  for 
the  effects  of  attenuation  due  to  pore  crushing  in  tuff.  This 
happens  because  the  high  frequency  spectra  are  dominated 
by  cavity  reverberations.  This  is  true  for  overdecoupled, 
decoupled,  and  partially  coupled  explosions. 
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